Phase-sensitive tests of the pairing state symmetry in Sr2Ru04 
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Exotic superconducting properties of Sr2Ru04 have provided strong support for an unconventional 
pairing symmetry. However, the extensive efforts over the past decade have not yet unambiguously 
resolved the controversy about the pairing symmetry in this material. While recent phase-sensitive 
experiments using flux modulation in Josephson junctions consisting of Sr2Ru04 and a conventional 
superconductor have been interpreted as conclusive evidence for a chiral spin-triplet pairing, we pro- 
pose here an alternative interpretation. We show that an overlooked chiral spin-singlet pairing is also 
compatible with the observed phase shifts in Josephson junctions and propose further experiments 
which would distinguish it from its spin-triplet counterpart. 



Unambiguous determination of the pairing state sym- 
metry is one of the key steps towards understanding the 
pairing mechanism in a continously growing class of un- 
conventional superconductors T|. Phase-sensitive exper- 
iments, capable of identifying angular dependence of the 
superconducting order parameter, have provided a cru- 
cial evidence for a dominant d-wave orbital symmetry in 
cuprate superconductors However, much less is 

known for other unconventional superconductors such as 
heavy fermions, charge transfer salts, and cobaltates. In 
particular, there is compelling evidence for an unconven- 
tional pairing in Sr2Ru04 |a|0|j with the strong possi- 
bility of spin-triplet superconductivity which would be a 
solid-state analog of superfluid He'^ 7] . 

In superconductors with inversion symmetry an order 
parameter (gap matrix) can be expressed as A(k) = 
Ao(k)ia'y for spin singlet and A(k) = a ■ d(k)z(Tj, for 
spin-triplet pairing. Here a are the Pauli spin matrices 
and scalar (vector) Aq (d) has even (odd) parity in the 
wavevector k. Often the symmetry of both the orbital 
and the spin part of A(k) remains to be identified and 
the lack of related understanding comes from the diffi- 
culty in performing phase-sensitive experiments. 

While numerous previous experiments probed the pair- 
ing symmetry of Sr2Ru04 [6j, in this context, re- 
cent phase-sensitive experiments that provide angle- 
resolved information are particularly important. The 
corresponding measurements were performed in a super- 
conducting quantum interference device (SQUID) geom- 
etry, consisting of a pair of Auo.5lno.5/Sr2Ru04 Joseph- 
son junctions. Since Auo.sIuq.s is a conventional s-wave 
superconductor, the observed modulation of critical cur- 
rent in an applied magnetic field was interpreted as con- 
clusive support for the phase shifts from an odd-parity 
spin-triplet pairing in Sr2Ru04 8, 9]. 

A similar SQUID geometry was initially proposed 
to study possible p-wave pairing in heavy fermions 
and later also used for identifying d-wave pairing in 
cuprates 0- The critical current Ic is modulated in the 
applied magnetic field as 

/,«cos($/$o + (512/2), (1) 

where $ is the flux threading the SQUID, $0 is the flux 



quantum, and 612 is the intrinsic phase shift of the order 
parameter between the two tunneling directions. For a 
conventional s-wave SQUID S12 = and Ic has a maxi- 
mum at $ = 0. In contrast, a phase shift S12 = tt, char- 
acteristic of unconventional pairing jl^ . yields a mini- 
mum /c at = 0. The modulation of external flux to- 
gether with the fabrication of junctions with varying tun- 
neling directions in SQUID geometry therefore provides 
an angle-resolved phase-sensitive information about the 
superconducting pairing symmetry. 

The suggested chiral n-wave (CpW) state with the 
triplet order parameter |l3|. 

d(k) oc {k^ + iky)z, (2) 

in which the spins of the Cooper pairs lie in the RUO2 
plane (_L d) , is indeed compatible with the experi- 
ment However, we show here that it is not the only 
candidate. There exists another pairing state, allowed by 
the tetragonal symmetry of Sr2Ru04, the singlet chiral 
d-wave (CdW) state ^Eg{c) with Ao(k) c>c {k^ + iky)kz, 
or, more accurately 

Ao(k) oc {kx + iky) sin fc^c, (3) 

which is equally consistent Ts'l with the phase shifts ob- 
served in Ref. j^. We use our findings to propose an ex- 
perimental test which would discriminate between CpW 
and CdW pairing symmetries. 

Could experimental and theoretical reasons be used to 
rule out CdW and favor only the CpW state? The two 
main arguments in favor of the CpW come from muon 
spin resonance and Knight shift experiments ^ J£J . The 
former indicate a time-reversal symmetry breaking be- 
low the transition temperature Tc, incompatible with the 
d2.2_j^2-wave state in cuprates, but fully compatible with 
either CpW or CdW symmetry. The Knight shift {K) 
was initially interpreted as firm evidence for a triplet 
state with in-plane spins (like CpW), since no change 
of the in-plane spin susceptibility below Tc was found. 

Even in singlet superconductors {e.g., vanadium), K 
could remain invariant below T^. Such behavior is usually 
attributed to (a) spin-orbit induced spin-flip scattering 
which suppresses the effect of the superconductivity on 
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K and (b) an accidental cancellation of the spin, dipole, 
and orbital contributions of the Fermi-level electrons to 
K which leaves only superconductivity-insensitive contri- 
butions such as the Van Vleck susceptibility. However, a 
quantitative analysis (171] shows that the spin-orbit cou- 
pling in Sr2Ru04 is too weak for scenario (a) while the 
accidential cancellation, required for scenario (b), does 
not occur Jl^]. Thus, neither of the two explanations of 
a constant K arising from singlet pairing is applicable. 

This would have made the Knight shift argument for 
CpW very convincing, if not for the recent experiment 
showi ng t he same result in a field perpendicular to the 
plane 'l^. It was proposed that the testing field 
of 0.02 T may be enough to induce a phase transition 
from the CpW in Eq. Q to a state with d||x. How- 
ever, this is highly unlikely: (i) the d||x state would have 
an additional horizontal line node, as compared to the 
dcx {kx + iky)z state and therefore lose a large part of the 
pairing energy ( ~ A per electron, A > 1.4 K ^ x0.02 
T); (ii) although in the d||x state the spins of the pairs lie 
in the yz plane, there is no y — z symmetry (as opposed to 
the xy plane) and it is not a priori clear whether the mag- 
netic susceptibility of the Cooper pair will be the same 
as for the normal electrons. Since d||x is not allowed for 
a tetragonal symmetry, it may only appear as a result 
of a second phase transition below Tc, which has never 
been observed in Sr2Ru04; (iii) the spin-orbit part of the 
pairing interaction, which keeps the spins in the xy plane, 
despite z being the easy magnetization axis _20j ■ would 
have to be weaker than ^bx0.02 T=1.1 ^eV=0.013 K, 
an energy scale much too small for the spin-orbit cou- 
pling in Sr2Ru04. So, neither the old theories for the 
lack of a Knight shift reduction below Tc, nor the new 
explanation in terms of the magnetic-field rotated order 
parameter withstand quantitative scrutiny; the Knight 
shift in Sr2Ru04 remains a challenge for theorists. Until 
this puzzle is resolved, we cannot use the Knight shift 
argument for the pairing symmetry determination. 

We now turn to the experiments of Ref. Q and com- 
pare Josephson tunneling between an s-wave supercon- 
ductor and cither (a) an even parity (spin singlet) su- 
perconductor or (b) an odd parity (triplet) superconduc- 
tor. In the first case, the Josephson current between 
a conventional s-wave superconductor and an unconven- 
tional spin-singlet superconductor, represented by the or- 
der parameters A^-wavc and Ao(k), respectively, can be 
expressed as [i^ 

Joc(TkIm[A:_Ao(k)])^^, (4) 

which depends on the relative phase between the super- 
conducting order parameters. The averaging is over all 
states at the Fermi surface (FS) where the Fermi velocity, 
vp, has a positive projection on the tunneling direction 
represented by the unit normal n (perpendicular to the 
interface plane, see Fig.^l and Tk is the tunneling prob- 
ability. For a thick rectangular barrier of width w and 



FIG. 1: Schematic sample and the Fermi surface geometry 
for phase-sensitive SQUID measurements from Ref. Q. (a) 
Auo.5lno.5/Sr2Ru04 junction geometry with an interface nor- 
mal n. (b) Possible deviation of n from the ab crystallographic 
plane, (c) Warping of the Sr2Ru04 Fermi surface. The mag- 
nitude of the Fermi wavevector k_F is generally different from 
the one corresponding to the cylindrical Fermi surface (kpo)- 
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height U |23| we can obtain 

16m^K^VLVRexp[—2Kw — k?,w/K] 



{k? + to^i;|^)(k^ -|- rri^vj^) 



(5) 



where n = ^2m{U — ji) such that wk ^ 1 (in the thick 
barrier limit), m is the free-electron mass, /i is the Fermi 
energy, and we set h = 1. We use vl,r to denote nor- 
mal components of the Fermi velocities in the two su- 
perconductors and fc|| is the component parallel to the 
interface. From Eq. Q we see that Tk is sharply peaked 
when vp \\ n. 

In the second case, the Josephson current between a 
singlet and a triplet superconductor becomes [21I Esf 



J cx (TkIm[A5(k)d(k) • (n x k)] 



FS ' 



(6) 



where we use Tk to denote that, unlike Tk, it con- 
tains matrix elements corresponding to spin-flip tunnel- 
ing, for example, due to magnetic interfaces or spin- 
orbit coupling. For nonmagnetic barriers and in the ab- 
sence of spin-orbit, there is no spin-flip scattering, there- 
fore Tk = and the Josephson current vanishes identi- 
cally Qmii. 

From Eqs. I@J and © we can directly infer that for 
c-axis tunneling, n || c, J = for both CdW and CpW 
states [J dkxdky(kx + iky) — 0]. For tunneling precisely 
in the ab plane, the current is also zero for CrfW while 
for CpW it only vanishes at n |j k. 

We consider a model of a quasi-two-dimensional (2D) 
layered superconductor which has a nearly cylindrical FS 
with a small c-axis dispersion originating from the weak 
inter-layer hopping. In Fig. ^a) we show the sample 
geometry used Ref. and in Fig. ^h) represent a warp- 
ing of the Fermi surface. In the first approximation for 
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Sr2Ru04 such a warping can be expressed as |2£ 



(^) 



2m 



[1 + £ COS kz 



(7) 



where |e| ^ 1 is the warping parameter (e « —7 x 10""* 
[2^), c is the lattice constant along the crystallographic 
^-direction, and ki?(z) is the z-dependent projection of 
the Fermi wavevector in the xy plane [see Fig.^c)]. It is 
convenient to resolve the Fermi wavevector in cylindrical 
coordinates {kp{z), if, kz) with 



kpiz) — fcFo/[l + £ cos kzc\ 



1/2 



(8) 



where kpo — (2to^)^/^ and for = [see Fig. ^c)] 
kpiz) kpx- 

Schematic representation of the SQUID geometry in 
Fig. ^a) (adapted from Ref. 0) is an oversimplifica- 
tion. While efforts were made to fabricate edges either 
precisely parallel or perpendicular to the c-axis, in the 
actual samples the direction of interface planes or their 
corresponding normals changes gradually from a to c- 
direction. In several samples j^j an interface nearly par- 
allel to the ab plane, at Auo.5lno.5/Sr2Ru04 junction, 
was covered by an insulating oxide [see Fig. QJa)]. It is 
then plausible to expect that the normal to such inter- 
face could deviate from the ab plane. In Fig. ^h) we 
depict a generalized situation in which an interface nor- 
mal, n = (np,n^,nz) with |nz| <g; 1, need not lie exactly 
in the crystallographic ab plane of Sr2Ru04. We show 
below that the analysis of phase-sensitive measurements 
in terms of the two small but finite parameters, e and Hz , 
can provide a qualitatively different interpretation from 
those which a priori assume £ = 71^; = 0. 

For a conventional superconductor with the FS larger 
than the one of Sr2Ru04, the Josephson tunneling across 
a thick rectangular barrier can be obtained from Eqs. Q 
and © as 



J cc dk(5(ek — M)vF-nexp(— fc|u;/2K) 

./ vj? ■n>0 



/ vj?-n>C 

X Im(fca; -|- iky) sin kzC, 



(9) 



where A;| — k|, — (kp ■ n)^, ^ mvl ^, and the projec- 
tion of the Fermi velocity in Sr2Ru04 along n is 

kpoll + e cos kzc]^^^ ki,r.ecamkzC 
v_F-n = n„ - - — ; — -n^. 



m 



2m [1 -I- s cos kzc] 



(10) 

For a thick barrier, the integration can be simplified by 
noting that the dominant contribution comes from /c[| = 
0. The right hand side of Eq. (jSJ, in the leading order in 
e and Uz, can be then reduced to y^iTK / wckpQnz{l — e), 
such that 



Jn ^ Anz{l - e). 



(11) 



where A characterizes the normal state barrier trans- 
parency. Thus, with a tilted interface {uz 7^ 0) there 



is a finite Josephson current even in the absence of any 
Fermi surface warping [e = 0). 

To verify that our findings of finite Josephson current 
in CdW state are not limited to the specific assumption 
of a thick rectangular barrier, we also consider the rather 
different case of a strong (5-function barrier. The corre- 
sponding transmission probability is [23 . |26| 



{VL + VR^ + 4C/2 ' 



(12) 



where v^^r are the normal components of the Fermi ve- 
locities in the two superconductors and U (^ v\ is 
the scattering strength. From Eqs. Q), (|12|l we obtain 

Joe dk5{e\^ — iJ.)vp-nlm{kx + iky) sinkzC, (13) 

J V F -nX) 

where, unlike in the case of a thick barrier, we perform (p 
and kz integration. In the leading order, the right hand 
side of Eq. (fT^ is —irkpf^enz, and yields 



Js 



'A euz 



(14) 



where again A characterizes the normal state trans- 
parency. In contrast to the thick-barrier result, the cur- 
rent now vanishes in the absence of FS warping. From 
Eqs. Hll|) . (|14|l one can conjecture that for a general case 
J « A{s - e), where < s < 1. 

The presence of small parameters e and Hz in 
Eqs. I|ll|) . (|14|) shows that the Josephson current in the 
CrfW state would be reduced as compared to the conven- 
tional SQUID with s-wave electrodes. However, the al- 
ternative picture, based on the CpW state, also contains 
small parameters which should be kept in mind when in- 
terpreting the experiment of Ref. . In addition to the 
small relative strength of the spin-orbit coupling (quan- 
tified by the admixture of Si into a nonrelativistic 6*1 
state, or the smn-orbit induced band shift relative to the 
band width |23), there can also be another small factor 
— a ratio of the lattice constant and the superconducting 
coherence length |^, approximately 6 x 10~^ |^. 

Results from Eq. (HTJ, (jHI) confirm that the CdW 
state could be compatible with the phase shifts observed 
in Ref. Furthermore, the azimuthal dependence of 
an order parameter coincides for both CdW and CpW 
states. While the proposed symmetry arguments 0, |^ 
exclude most of superconducting states allowed in the 
tetragonal symmetry |29l | , these arguments alone are not 
sufficient to unambiguously identify the odd-pairing of 
CpW state. Instead, to confirm that a CpW state has 
indeed been observed, one would need to accurately cal- 
culate the expected magnitudes of the Josephson current 
for both chiral states. In particular we propose a mod- 
ification of the experimental configuration [Sj such that 
the interface plane would be slanted at w 45° with the c- 
axis. If the corresponding ratio of the Josephson current 
to the normal state conductivity becomes substantially 



4 



larger {riz is no longer small) than in Ref. [Sj, it would 
be strong support for the chiral singlet state in Eq. Q. 

Another important distinction between the two chi- 
ral states is the presence of nodes in the superconduct- 
ing gap. In contrast to the CpW state, CdW requires 
by symmetry a horizontal line node [see Eqs. l(2Jl, ®]- 
The idea of a horizontal line node has been enter- 
tained by experimentalists |211| and theorists [sJl for a 
while, although recently it has fallen out of favor. Still, 
some researchers insist on the existence of a horizontal 
line node |32l |. Moreover, in the Josephson experiments 
of Ref. evidence was found for a substantial de- 
pendence, albeit not necessarily for horizontal nodes, of 
the order parameter in Sr2Ru04 jl^. How could such 
a material with nearly 2D electronic structure develop a 
highly 3D superconducting state? To answer this ques- 
tion we point out the following facts: (a) practically no 
ferromagnetic spin fluctuations, favorable for a p-wave 
pairing, have been experimentally found in Sr2Ru04; (b) 
antiferromagnetic spin fluctuations at q= (2/3, 2/3, g^) 
have negligible z dispersion; (c) the crystal structure of 
Sr2Ru04, as opposed to its electronic structure, is fairly 
3D, so one can expect the electron-phonon coupling to be 
quite 3D as well (d) there is a sizeable O isotope effect 
in Sr2Ru04, which strongly changes with introduction 
of pair-breaking symmetries 34] . While electron-phonon 
coupling per se can only induce an s-wave pairing, such 
a pairing would be prevented by the strong antiferromag- 
netic spin fluctuations. However, for the proposed CrfW 
state, any 2D interaction cancels out, including the mag- 
netic interaction. Should the electron-phonon coupling 
have a maximum say, at (1/2,1/2,1/2), as opposed to 
(1/2,1/2,0), the CrfW state would have been immedi- 
ately stabilized providing a plausible scenario for spin- 
singlet superconductivity in Sr2Ru04. 

In conclusion, we have revealed that a completely over- 
looked chiral d-wave pairing state in Sr2Ru04 is equally 
compatible with the existing body of experimental data 
as the generally accepted chiral p-wave state. We have 
proposed phase-sensitive experiments in a SQUID geome- 
try with a variable tilting angle capable of unambiguously 
distinguishing between the two chiral states. 
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